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The Mars Atmosphere and Volatile EvolutioN (MAVEN) mission will 
determine the history of the loss of volatiles from the Martian 
atmosphere from a highly inclined elliptical orbit. MAVEN will launch 
from Cape Canaveral Air Force Station on an Atlas-V 401 during an 
extended 36-day launch period opening November 18, 2013. The 
MAVEN Navigation and Mission Design team performed a Monte Carlo 
analysis of the Type-II transfer to characterize; dispersions of the arrival 
B-Plane, trajectory correction maneuvers (TCMs), and the probability of 
Mars impact. This paper presents detailed analysis of critical MOI event 
coverage, maneuver constraints, AV-99 budgets, and Planetary 
Protection requirements. 


INTRODUCTION 

The Mars Atmosphere and Volatile EvolutioN (MAVEN) mission will determine the history 
of die loss of volatiles from the Martian atmosphere to space and in so doing will determine the 
impact of this process on the climate, the habitability, and the presence of liquid water. 1,2,3 
MAVEN will achieve this by providing a comprehensive picture of the present state of the upper 
atmosphere and ionosphere of Mars and the processes controlling these regions from a highly 
inclined elliptical orbit. 

MAVEN will launch from Cape Canaveral Air Force Station (CCAS) on an Atlas-V 401 
during a unique extended 36-day launch period which opens November 18, 2013. The worst-case 
combination of injection energy (C3) and outbound declination provides a substantial mass 
margin for the launch capability. Upon arrival at Mars on September 22 nd , 2014 (September 28 th , 
for a nominal launch window close), MAVEN will perform a designed pitched-over capture 
maneuver with a nominal thrust of 1321N, imparting a delta-V (AV) of approximately 1233 
m/sec over 36 minutes to achieve a 35-hr capture orbit at 75° inclination with a post-capture 
periapsis radius of 3779 km. Subsequent transition maneuvers will reduce the period to 4.5-hrs 
with the periapsis altitude within a required atmospheric density corridor of 0.05 kg/km 3 to 0.15 
kg/km 3 . A northern approach with a 5:45 pm local arrival time provides Earth visibility during the 
insertion bum. 
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As part of the type-II transfer trajectory design, the Mars planetary protection requirements 
must be satisfied for any uncorrected injection error over the entire launch period. Using United 
Launch Alliance (ULA) injection correlation data produced using ini tial GSFG and JPL tr ansfe r 
trajectory designs, the MAVEN Navigation and Mission Design team performed a Monte Carlo 
analysis to characterize; dispersions of the arrival B-Plane, the trajectory correction maneuvers 
(TCMs) which eliminate these dispersions, and the probability of Mars impact to assess planetary 
protection considerations. The optimized TCM-1 will mitigate operational impacts but is 
constrained by a minimum 4-m/s AV magnitude at the first trajectory correction opportunity. This 
first correction maneuver is used to exercise the Mars Orbit Insertion (MOI) thruster assembly. 
Based on analysis and the maneuver constraint, a single B-Plane aim-point was chosen for the 
entire 36-day launch period. 

Mission Overview 

The MAVEN mission will be the first mission devoted to understanding the Martian upper 
atmosphere and addressing the compelling questions regarding its initial formation and 
subsequent evolution. 4,5,6 These questions address the nature and history of Martian habitability 
by microbes, and how and why it has changed through time. As such, these questions fit cleanly 
into the Mars exploration program, whose broad goals include understanding the history of 
habitability and whether any organisms have ever existed on the planet. The MAVEN mission 
will explore the upper atmosphere of Mars and help scientists determine what role that loss of the 
atmosphere to space played in the history of the Martian atmosphere and climate. Over the last 
decade, data collected from many successful Mars missions support the view that liquid water 
existed at the surface early in Martian history. Increasingly, evidence points to loss of gases out 
the top of the atmosphere to space as an important, and possibly the dominant, process in the 
changing climate. Although recent measurements provide compelling evidence that loss to space 
has occurred, they do not allow a unique estimate to be made of how much gas has been lost, or 
to determine the specific processes by which the loss occurred. 


MAVEN measurements can be considered 
from two different perspectives. From the 
perspective of the “science goals”, it will make 
three different types of measurements. First, it 
will determine the present-day composition and 
structure of the upper atmosphere. Second, it 
will determine the present-day escape rate of gas 
from the upper atmosphere to space. Third, it will 
make measurements that allow us to extrapolate 
this escape rate to past times, when the solar 
wind and the solar ultraviolet fight (that drive the 
escape) were greater, to estimate the total amount 
of gas that has been lost. 

From an “observational” perspective, MAVEN takes three types of measurements. 5 ' 6 First, it 
will measure the properties of the upper atmosphere as the spacecraft passes through the upper 
atmosphere on each orbit. These allow a very detailed look at one place in the atmosphere on 
each orbit, and allow determination of the basic state of the upper atmosphere. Second, it will 
make remote-sensing measurements of a large part of the planet from the high-altitude parts of its 
orbit. This will allow the point measurements to be extrapolated to global conditions, and will 
provide a good understanding of the geographical variations that can take place. Third, MAVEN 
will measure the energy inputs into the upper atmosphere that drive the processes that lead to 
escape. This will include the properties of the solar wind as it hits Mars, of solar ultraviolet fight, 
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and of solar storms, all of which can affect the behavior of the top of the atmosphere. Figure 1 
shows the different atmospheric loss and energy processes that MAVEN will measure. Neutral 
processes are shown in blue, ion and plasma processes in red, and solar energetic inputs are 
shown in the upper left. 

Spacecraft Overview 

The MAVEN spacecraft, shown in Figure 2, is the latest in a series of Lockheed Martin Mars 
orbiters to be developed for NASA (past orbiters include the Mars Global Surveyor, Mars 
Odyssey, and Mars Reconnaissance Orbiter (MRO)). 5 Spacecraft features include redundant star 
trackers and Inertial Measurement Units (IMU); two 50-Iiour nickel hydrogen batteries; two fixed 
solar arrays; low-gain and high-gain antennas; a mono-propellant system; redundant Command 
& Data Handling system; and a fault-tolerant Mars Orbit Insertion (MOI) design. In addition, it 
is a 3-axis stabilized sun-pointing spacecraft, with a fixed high gain antenna. The spacecraft 
accommodates eight instruments (see Figure 3) that are segregated into three distinct packages: a 
Particles & Fields (P&F) package being developed by UC-Berkeley/Space Sciences Laboratory; a 
Remote Sensing Package being developed by CU-LASP; and the Neutral Gas and Ion Mass 
Spectrometer (NGIMS) being developed by GSFC. 
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Figure 2. MAVEN Spacecraft 

MAVEN TRANSFER TRAJECTORY DESIGN 

MAVEN will launch from Cape Canaveral Air Force Station (CCAS) during a 36-day launch 
period (a 20-day nominal period with 16-days extended) opening November 18, 2013 on an 
ATLAS-V 401. The worst-case combination of injection energy and launch declination results in 
an allowable launch mass margin for the proposed 2354 kg current best estimate wet mass. Upon 
arrival at Mars on September 22, 2014 (September 28, for a Dec 7, 2013 nominal launch period 
close), MAVEN will insert into a capture orbit with a 35-hour period. Figure 3 shows the transfer 
trajectory and the dates of launch and arrival at Mars including the locations of all the planned 
TCMs. As with most planetary launches, there is a trade between the launch energy (C3), launch 
declination, and the arrival V-Infinity as these parameters vary over the launch period and they 
can become the driver for the mission design and propulsion system selection. 
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Figure 3 MAVEN Type - II Ballistic Transfer 


For MAVEN, a change in arrival date by a 
week can result in an increase in the arrival V- 
infinity without a significant change to C3. 
The increased V- infinity yields an increase in 
MOI AV of ~ 50 m/s. Since the propellant 
tank is full, a late arrival scenario is highly 
unfavorable as there is no additional 
propellant margin. Figure 4 shows contours of 
the C3 and V-Infinity for the desired launch 
period. One can see in Figure 4a that the C3 
over the launch departure dates are relative 
flat and decrease over the period by 2.6 
km 2 /s 2 . The V-Infinity in Figure 4b however, 
has only a short duration where it re mains flat 
near the minimum of 3.17 km/s. With a 
decreasing C3 over the launch period, the 


MOI AV determined from V-infinity determines the close of the launch period. The launch 
declination remains under 28 degrees over the nominal launch period and thus does not pose a 
concern for reduced mass capability. Table 1 provides launch period information for the MAVEN 
mission. 
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Figure 4 (a) and (b). C3 and C3 / V-Infinity Contour Plots 


Table 1: Analytical Launch Period Parameters 


Launch Period 

Open 

Mid 

Close 

C3 (km 2 /s 2 ) 

1207 

10 30 

9.40 

Declination (deg) 

13.38 

18.84 

26.93 

Right Ascension (deg) 

198.26 

200.58 

200.88 

Arrival date 

9/22/2014 

9/24/2014 

9/28/2014 

Arrival V-Infinity 

3 17 

3 15 

3.17 
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Selection of B-Plane Components 

The selection of the target B-Plane begins with the completion of the determination of the 
ballistic type-II transfer. We initiated simulations using a full ephemeris and high fidelity 
environmental models to define transfer trajectories that propagated to the Mars periapsis. The B- 
plane components as defined in Figure 5 are computed at the epoch of the periapsis event. B- 
plane components are defined similarly to past Mars mission design references. The coordinates 
used for the B-Planes in this paper are based on the Mars Mean Equator IAU definition, with k 
being the Mars’ North Pole. Note that the B.R negative direction is ‘upward’ and the B-angle 
increases in a counterclockwise direction. 

The selection of the targeted B-Plane components is a result of several considerations. It must 
meet the planetary protection goals of less than le-4 probability [reference 9] and meet the 
operational goal of a first TCM of 4 m/s minimum. This minimum AV is chosen so that the 
designed propulsion system can reach a steady state condition to successfully test the orbit 
insertion thrusters at an early portion in the transfer timeline. The current schedule has the first 
maneuver, TCM-1, occurring 15-days after injection. The chosen schedule permits time for 
navigation and other spacecraft operations to be completed. 



Transfer Injection Position State Generation 

To obtain the initial states for our process, a simplified B-plane determination process was 
used and iterated with KSC and the Launch Vehicle Provider, ULA. Using high fidelity 
simulations with environmental models from all past Mars missions and estimates of the initial 
launch conditions, we targeted to the desired orbit insertion conditions. Starting with this 
information and an esti m ated TIP error matrix based on preliminary information exchanged with 
ULA, the direction of the maximum and minimum TCM-1 AV vectors can be computed for a 
preliminary estimate of the target B-Plane components that meets the planetary protection and 
TCM-1 AV requirement. From this iterated analysis, a selected B-Plane aimpoint was chosen 
with components of B.T =16100 km, and B.R = -19300 km. 

This aimpoint applies to the optimal condition of the daily launch window for the full launch 
period. For the first and second Confi g uration, Performance and Weight Status Report (CPWSR) 
delivery the B-Plane targets were held constant across the nominal launch period. Initial targeting 
simulations gave the results shown in Table 2 over the entire launch period. The table presents the 
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launch energy (C3) and outgoing asymptote angles of DLA and RLA for each date. Figure 6 
shows the trending and the variation in C3, DLA, and RLA associated with discrete arrival times. 


Table 2 — MAVEN Launch Ini 


urination for Nominal and Extended Period 


Launch 

Opportunity 

Launch Date 

IIP Epoch 
(LTCG) 

C3 

(km A 2/sec A 2) 

DIA 

(deg) 

HI -A 
(deg) 


Launch 

Opportunity 

Launch 

Date 

Epoch 

(LICG) 

C3 

(km A 2/sec A 2) 

DL4 

(deg) 

HI A 
(deg) 

1 

11/18/13 

20.38.58 

12 23 

12.78 

198 26 


21 

12/08' 13 

16.55.1 J 

9 31 

26.72 

200.35 

? 

1M9/13 

20.32.25 

1198 

i3.24 

198 52 


22 

12/09/13 

16 39-48 

9 30 

27 42 

200,87 

3 

11/20/13 

20:25.33 

11 73 

13, .'2 

198 78 


23 

12/10/13 

16.04:21 

9 35 

28 50 

200.84 

4 

1 1/2 1 ' 1 3 

20 18:22 

il 50 

14.24 

199 03 


24 

12 1 M3 

15.23 36 

9.41 

29.60 

200.83 

5 

1 1 '22/1 3 

20,10:48 

1127 

14 80 

199 27 


25 

12/12-13 

J7 29 38 

9 53 

30.88 

200.85 

6 

1 1/23/13 

20:04:30 

11 08 

15.24 

199.61 


26 

12' 13/13 

17,04-34 

9.66 

32 07 

200.89 

7 

11/24/13 

19:57 54 

10.89 

15 70 

199 94 


27 

12,14/13 

16 49:49 

9 73 

32.83 

200.93 

8 

11/25/13 

19-48 51 

10.69 

lo 39 

200.16 


28 

12/15 13 

16 3/08 

9.91 

34.09 

201 Ul 

9 

11/26 13 

19 38:42 

10.1-9 

17 17 

200 33 


29 

12/16/13 

16 32.00 

10.02 

34 87 

201.07 

10 

11/27/13 

19:32.03 

1033 

17 57 

200 58 


30 

12/17/13 

16 46 Id 

10.i4 

35.67 

201.14 

11 

11/ 28/13 

19 21.30 

10 15 

18 28 

200.57 


31 

12*18/13 

16 43 00 

10,27 

36.46 

201,21 

12 

11 29/13 

19 1 1.02 

9 98 

18 98 

200.63 


32 

12/19/13 

16 I® 06 

10 55 

37,82 

201 35 

13 

11.30/13 

18:5948 

9.84 

19 ^5 

200 70 


33 

12,20 13 

16.30.25 

10^3 

38.64 

201.43 

14 

12/01/13 

18.47 36 

9 n l 

20 57 

200 76 


14 

12/21/13 

16 36.01 

11 10 

4010 

20 1 62 

15 

12/02^13 

18.27.33 

9 60 

21 18 

200 39 


35 

12/2213 

16:26 48 

1L05 

40.04 

201 57 

16 

12/03/13 

18 23 27 

9 50 

22 08 

200 90 


36 

12/23/ 13 

16:25.50 

10.91 

39.69 

201 38 

17 

12/04/13 

18 13.23 

9 42 

23.02 

200.90 







18 

12/05/13 

17:56.42 

9 37 

23 99 

200 38 







19 

12WI3 

17:38 13 

9 33 

24.99 

200 85 







20 

12/07/13 

1723 17 

9 30 

25 69 

200.83 
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Figure 6. MAVEN C3, DLA and RLA over the Nominal and Extended Launch Period 


Covariance Analysis Process 

Covariance analysis of the type-II ballistic trajectory determines both the TCM AV budget 
and provides information on meeting the planetary protection requirements. Analysis began with 
the provided ULA Trajectory Injection Position (TIP) states. These states represent the post 
separation conditions without errors from the centaur injection. A set of figure of merit matrices 
was also provided with each state as part of the ULA CPWSR-1 delivery. We sampled a 
probability distribution associated with each eigenvalue and eigenvector by using the following 
procedure. 
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Monte Carlo analysis process 

• Ingest correlation & covariance m atrix in the RTN frame 

• Multiply matrix for unit conversion (ft/sec to m/s) and copy covariance data to lower 
triangle 

• Use initial simulation state (TIP state) to adjust Covariance matrix into EME2000 

• Compute Eigenvalues and Eigenvectors using Matlab® ‘eig’ function 

o [V,D] = eig (C) where C is the 6x6 error covariance 

• Compute standard deviation (a) using square root (sqrt) of eigenvalues. We assume that 
perturbations in each direction follow a zero-mean Gaussian distribution. The standard 
deviation (a) of each distribution is. 

o std = sqrt (diag of D) 

• Use Matlabs’s normally distributed random number generator with the standard deviation 
values to scale the random numbers for our application. We define a realization vector 
(5), which is a change in the nominal injection state in Cartesian elements. This 
realization vector is found using; 

o S=V*T 

« Where V is the 6x6 matrix and T is a 6x1 matrix, of the random numbers generated using 
the standard deviation values defined in (a). This realization is added to the nominal 
Cartesian injection state and the injection epoch is defined as a constant value given by 
the ULA. That is to say, the off-nominal states are computed by addition of random 
errors to baseline state such that the new State = old state + error 

[Note: In Matlab this is scripted as dS(:,i) = randn(l,l)*Std(l)*V(:,l) + 

randn( 1 ,l)*Std(2)*V(: ,2) + randn(l,l)*Std(3)*V(:,3) + randn(l,l)*Std(4)*V(:,4) + 

randn(l,l)*Std(5)*V(:,5) + randn(l,l)*Std(6)*V(:,6);] 

• State is loaded into AGI’s Astrogator® and used in propagation and/or targeting for 1000 
samples 


Covariance Analysis Results 


The overall design process is 
to select the initial guess for a 
biased aimpoint (B.T, B.R) and 
determine die dispersion at Mars | 
arrival epoch. In doing so, we * 
produced an injection covariance 
to map the encounter analytically 
and generate a B-Plane covariance 
using Monte Carlo runs. We use a 
Probability Density Function 
(PDF) of a calculated distribution t 
(assumed linear with aim-point) to j 
determine the P[Impact] for a 
given aimpoint. From this we can 
then determine the AV required to 
correct a given aimpoint bias. The 
challenge for MAVEN is to find a 
B.T and B.R set that minimizes AV. 



Figure 7. Sample Histogram of Post Injection Errors 

constrained such that AV > 4 m/s and P[I] < le-4. These 
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results were confirmed by independent verification between AGI’s Astrogator Module and JPL’s 
MONTE software as part of the standard MAVEN navigation support practice. Figure 7 shows a 
sample histogram of the generation of the errors on the TIP state using the above process. The 
figure shows the position and velocity errors in the EME2000 frame. Note that the 3-sigma errors 
in position are approximately 10 meters while the velocity errors are 15 m/s. 
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Figure 8 Sample B-Plane Distributions for Optimal Nov 28 th 2013 Launch 

Figure 8 presents a sample B-plane distribution for the MAVEN trajectory if launched at the 
opening. As seen the angle of the distribution for this launch opening case is near 12-deg and is a 
function of the launch date and the orientation of the relative orbit planes of Earth and Mar s 
From the initial states, with errors, the distribution of the B-plane components can be used to 
determine the PDF and therefore the probability of impact with Mars. 



Probability of Mars Impact 

Equations 1-4 provide the background of the method used to compute the impact probability. 
Equation 1 shows the calculation of the 2-D Gaussian PDF. Equation 2 shows the probability of 
being in area [A], and equation 3 shows the calculation of the ‘B-Plane’ area that would impact 
the Martian surface. Thus, impact probability can be calculated using equation 4. 




(i) 


'M-lf/W* 

A 


( 2 ) 


T- 

\mp 


= r ji + 2-GM km 

1 C 3 r„ 


(3) 


P{I\ = ]dr]d6(rf{r,e)) 
0 0 


(4) 


Based on these equations and the data from the high fidelity simulations using AGI’s 
Astrogator 18 '' (V9.2) and Matlab®, the data shown below is generated. Figure 9 shows a sample of 
the probability of impact with the highlighted contours (blue) indicating the le-4 requirement. 
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Also shown is the 4 m/s TCM-1 AV 
magnitude (red) elliptical shaped contour 
to meet the insertion orbit requirements. 
Several choices can be made regarding 
the selection of the aimpoint that meets 
both requirements. Based on these data, 
impact and TCM-1 AV cost can be 
computed and compared across the 
launch period or launch day. 

Launch Period B-Plane 
Distribution and Impact 
Probability 


Following the procedures and 
equations shown above, Monte Carlo „ 

analysis was completed for the launch Fl S ure 9 * Combined Probab,llt y of Impact 

period opening, mid, close and an and TCM-1 AV (m/s) Magnitude 

extended launch period date. This analysis represents the launch dispersion estimates for the 
CPWSR-1 delivery and is currently undergoing an update. As the analysis was completed it was 
recognized that the angle of the B-Plane distribution ellipses’ semi-major axis with respect to the 
B.T direction (in the Mars orbit plane by definition), increased both as the launch date moved 
from launch opening to launch closure and also minimally across the daily 2-hour launch 
window. The major and minor axes amplitudes also varied over the period with an average of 2e5 
km and 2e4 km respectively. The distribution geometry is a function of the angle between the 
Earth and Mars orbit planes, launch geometry, and launch injection errors. 



The next series of figures presents the B-plane distribution and related probability contours 
for launch open (Nov 18*), mid (Nov 27*), close (Dec 07*), and extended (Dec 15 th ) of the 
launch period. Also included is the effect of a 2-hr launch window for each launch date, as shown 
in the figures from left to right. For epochs between those shown, we found a linear effect. 
Figures 10 through 13 shows the B-Plane distribution for these launch dates while Figure 14 
through 16 shows the probability of impact. Information for an extended launch date of Dec 15* 
is provided to show the continuation of the analysis for evaluation of launch dates past the 
nominal period end. 


The Probability of Impact, P[I], for the shown B-Plane distributions all meet the required 
P[I] < le-4. The way to interpret the probability plot contours is similar to the B-Plane plots. A 
small circle is shown that represents the nominal aimpoint which meets the P[I] and TCM-1 AV 
requirements and also represents the targeted point for all ULA deliveries. The P[I] contour lines 
represent the probability of impact of each possible B-Plane aimpoint on the figure based on 
equations 1-4 with the mean being the center of the distribution and r imp being the impact radius. 
Figure 17 shows the probability of impact over the entire launch period and that the least 
favorable P[I] occurs at the launch opening. 
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Figure 10. B-Plane Distribution for Launch Open, Nov 18 2013, with 2-hr window 
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Figure 11. B-Plane Distribution for Launch Mid, Nov 27 th 2013, with 2-hr window 
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Figure 12. B-Plane Distribution for Launch Close, Dec 07 2013, with 2-hr window 
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Figure 13. B-Plane Distribution for Launch Extended, Dec 15 th 2013, with 2-hr window 
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Figure 15. B-Impact Probability for Launch Mid, Nov 27* 2013, with 2-hr window 
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Figure 15* B-lmpact Probability for Launch Close, Dec Or 2013, with 2-hr window 
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Figure 16. B-Impact Probability for Launch Extended, Dec 15 th 2013, with 2-hr window 
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Figurel7. Impact Probability over Launch Period 





TCM-1 AY Statistics 


With the B-Plane distribution statistical analysis completed, we determined the TCM-1 
budget for the transfer trajectory. For each of the cases provided above, the distribution of the 
TCM-1 AV performed at injection plus 15 days is shown. As can be seen the variation in the 
TCM-1 AV magnitude can be quite large, ranging from approximately 2 m/s to over 15 m/s 
depending on the launch date and the time of day of launch. Figures 18 through 21 show the 
TCM-1 AV magnitude (units of km/s) distribution for each launch condition corresponding to the 
above B-plane distribution and P[I]. Using the AV magnitude and MATLAB statistical functions 
(e.g. percentile, mean, median). Table 3 provides the statistical analysis of these plots and the 3- 
sigma and AV-99 for the nominal launch period. The extended launch period date of Dec 15* 
indicates a AV magnitude that approaches 50 m/s. While the AV magnitude is larger than the 
allowable TCM-1 AV budget of 20 m/s, it is still within family when considering the total AV 
budget and the desire to launch during the 2013 cycle. The epoch of TCM-1 can also be adjusted 
to vary the magnitude of the TCM-1 AV in order to maintain a minimum 4 m/s requirement. 



Figure 18. TCM-1 AV s for Launch Open, Nov 18* 2013, with 2-hr window 
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Figure 19. TCM-1 AV s for Launch Mid, Nov 27 th 2013, with 2-hr window 



TCM-1 Dtf magiMude* 15 <&ys TCM-1 CM magntude at IS daps TCM-1 OT magnikjOe at 15durs 


Figure 20 TCM-1 AV s for Launch Close, Dec 07 th 2013, with 2-hr window 



Figure 21 TCM-1 AV s for Launch Extended, Dec 15 th 2013, with 2-hr window 
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Table 3. TCM-1 Statistical AV In 


ormation 


Launch date and window 

Max 

(m/s) 

Mirt 

(m/s) 

Mean 

(m/s) 

Std 

(m/s) 

3-sigma 

(m/s) 

DV-99 

(m/s) 

Nov 18 Minus lhr 

9.6 

1.2 

4.9 

1.2 

8.6 

7.9 

Nov 18 Optimal 

12.8 

1.7 

5.5 

1.5 

10.1 

9.6 

Nov 18 Plus lhr 

12.9 

2.4 

6.6 

1.8 

^11.9 

11.2 

Nov 27 Minus lhr 

9.9 

2.1 

5.4 

1.2 

9.0 

8.6 

Nov 27 Optimal 

13.0 

2.3 

6.4 

1.7 

11.4 

10.5 

Nov 27 Plus lhr 

15.1 

2.7 

8.2 

2.0 

14.1 

13.0 

Dec 07 Minus lhr 

11.1 

1,8 

5.7 

1.5 

10.3 

9.7 

Dec 07 Optimal 

15.5 

3.8 

8.9 

1.9 

14.6 

13.4 

Dec 07 Plus lhr 

20.3 

6.8 

13.1 

2.3 

20.0 

18.6 


ADDITIONAL DESIGN CONSIDERATIONS 
Modeling and Optimization 

The environmental models used for the analysis presented above are representative of all 
previous Mars missions and the MAVEN Mission Design and Navigation team has thoroughly 
verified consistent and accurate trajectory designs using both the AGI Astrogator V9.2 and JPL’s 
MONTE software. As part of the ongoing analysis the team is well aware that optimizing a 
combination of TCM-1 AV at 15-days after launch with the TCM-2 AV performed 90 days after 
launch would benefit the overall TCM budget as well as permit more operational options. The 
TCM-1 in this optimization process could be held to the required 4 m/s magnitude. 

Orbit Plan Angle 

As seen in several figures, the slope (angle wrt the B.T direction) of the B-Plane distribution 
ranges from ~1 1-deg to ~45-deg over the launch period. This angle plays an important function in 
the magnitude of the TCM-1 AVs. The importance of the distribution is in the TCM AVs to 
correct the trajectory from the biased aimpoint back to the required insertion conditions. The AVs 
are minimal in the B.T direction as that direction is a heliocentric orbit period correction. The AV 
to change the B.R direction is larger as it is associated with a heliocentric inclination or arrival 
declination change. Also the location of the TCM-1 maneuver is not optimal wrt to inclination 
changes nor represents an optimal broken plane maneuver. With the opening launch date arrival 
conditions, the heliocentric inclination is ~ 2-deg and increases nearly 1-deg towards the launch 
closure date. The Mars incoming declination varies from -20-deg at opening to -33-deg at closure. 

MOI constraints 

As part of the various rules and goals with respect to observation of critical events, the timing 
.of the MOI and the targeted Mars arrival epoch was controlled such that two independent DSN 
sites could maintain coverage of the complete MOI duration. This is a straight forward targeting 
strategy in that the arrival epoch of Mars periapsis is targeted in the overall design process. This 
arrival time is dependent upon the dual coverage of both Canberra and Goldstone DSN sites for 
MOI and is centered at 02:00 hrs. Since the MOI maneuver is approximately 36 minutes in 
duration, even with fault tolerant measures in place, this dual coverage permits either site 
coverage throughout the maneuver. Other launch period dates are similar since file arrival time is 
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fixed and any launch injection or navigation errors can be fixed with s m all correction maneuvers 
of approximately 10 cm/s in magnitude. 

CONCLUSION 

We have described the methodology that was used to perform the MAVEN injection 
dispersion analysis. This methodology includes running a Monte Carlo analysis based on a single 
injection opportunity that can be assumed to be independent of other errors and is considered the 
baseline mission design case. From this analysis we confirm that the MAVEN project budgets 20 
m/s for the TCM-1 AV that is performed at 15-days after injection for 3a injection dispersions. 
This recommendation is equivalent to the 3a of the TCM-1 magnitude AV with a 50% margin at 
launch period opening and 25% margin at launch period close. Changes to the MOI maneuver are 
negligible with respect to the results of the Monte Carlo analysis. The TCM-1 AV magnitudes 
depend upon the specific Atlas-V 401 covariance matrix delivered. The analysis presented will be 
repeated using ULA CPWSR-2 and CPWSR-3 data for a definitive MAVEN AV budget. 
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